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ABSTRACT 

We present computer simulations of interferometric imaging with the Very Large Telescope Interferometer (VLTI) of 
the European Southern Observatory (ESO) and the Astronomical MultiBEam Recombiner (AMBER) phase-closure 
instrument. These simulations include both the astrophysical modelling of a stellar object by radiative transfer 
calculations and the simulation of light propagation from the object to the detector (through atmosphere, telescopes, 
and the AMBER instrument), simulation of photon noise and detector read-out noise, and finally data processing 
of the interferograms. The results show the dependence of the visibility error bars on the following observational 
' parameters: different seeing during the observation of object and reference star (Fried parameters tq, object and ro^rof. 
[t I I ranging between 0.9 m and 1.2 m), different residual tip-tilt error ((5tt, object and (5tt,ref. ranging between 0.1% and 20% 
of the Airy disk diameter), and object brightness (_ft'object=0.7mag to 10.2mag, ii'ref.=0.7mag). Exemplarily, we 
focus on stars in late stages of stellar evolution and study one of its key objects, the dusty supergiant IRC -1-10 420 
T-H ■ that is rapidly evolving on human timescales. We show computer simulations of VLT interferometry (visibility and 
' phase closure measurements) of IRC -1-10420 with two and three Auxiliary Telescopes (ATs; AMBER wide-field 
. mode, i.e. without fiber optics spatial filters) and discuss whether the visibility accuracy is sufficient to distinguish 
' between different theoretical model predictions. 
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O ■ 1. INTRODUCTION 

''^ ' The Very Large Telescope Interferometer^ (VLTI) of the European Southern Observatory (ESO) with its four 8.2 m 

I ' unit telescopes (UTs) and three 1.8 m auxiliary telescopes (ATs) will certainly establish a new era of optical and 

O infrared interferometric imaging within the next few years. With a maximum baseline of up to more than 200 m, the 
VLTI will allow the study of astrophysical key objects with unprecendented resolution opening up new vistas to a 

^ [ better understanding of their physics. 

^ I The near-infrared focal plane instrument of the VLTI, the Astronomical MultiBEam Recombinerii (AMBER), 

• ""j ' will operate between 1 and 2.5 fim and for up to three beams allowing the measurement of closure phases. In a second 

rS . stage its wavelength coverage is planned to be extended to 0.6 /im. Objects as faint as K — 20 mag are expected to 

■ be observable with AMBER when a bright reference star is available, and as faint as K — 14 mag otherwise. 

Among the astrophysical key issued are, for instance, young stellar objects, active galactic nuclei and stars in 
late stages of stellar evolution. The simulation of interferometric imaging of a stellar object consists in principle of 
two components: 

(i) the calculation of an astrophysical model of the object, typically based on radiative transfer calculations 
predicting, e.g., its intensity distribution. To obtain a robust and non-ambiguous model, it is of particular 
importance to take diverse observational constraints into account, for instance the spectral energy distribution 
and visibilities. 

(ii) the determination of the interferometer's response to this intensity signal, i.e. the simulation of light propagation 
in the atmosphere and the interferometer. 

Often, only one of the above parts is considered in full detail. The aim of this study is to combine both efforts and to 
present a computer simulation of the VLTI performance for observations of one object class, the dusty supergiants. 
For this purpose, we calculated a detailed radiative transfer model for one of its most outstanding representatives, 
the supergiant IRC -1-10 420, and carried out computer simulations of VLTI visibility measurements. The goal is to 
estimate how accurate visibilities can be measured with the VLTI in this particular but not untypical case, to discuss 
if the accuracy is sufficient to distinguish between different theoretical model predictions and to study on what the 
accuracy is dependent. 
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2. THE SUPERGIANT IRC +10 420: EVOLUTION ON HUMAN TIMESCALES 

The star IRC +10420 (= V|1302 Aql) is an OjUtstanding object for the study of stellar evolution. Its spectral type 
changed from FSI^ in 19733 to mid-A todayO'D corresponding to an increase of its effective temperature of 1000- 
2000 K within only 25 yr. It is heavily obscured by circujnstellar dust due to strong mass loss with rates typically 
of the order of several lO~''M0/yr {Mq: solar mass).lal3 IRC +10420 is believed to be a massive luminous star 
of initially ~ 20 to dCLMi^ currently being observed in its rapid transition from the red supergiant stage to the 
Wolf-Rayet ddphase.BE3EJi3'Ll Wolf-Rayet stars, in turn, finally evolve into a supernova explosion. IRC +10 420 is 
the only object observed until now in its transition phase to the Wolf-Rayet stage. 

Several infrared speckle and coronographic observati£jng|li0lll00 were conducted to study the dust shell of 
IRC +10420. The most recent study of Blocker et al.E3 reports the first diffraction-limited 73 mas (mas: milli- 
arcsecond) bispectrum speckle interferometry of IRC +10 420 and presents the first radiative transfer calculations 
that model both the spectral energy distribution and the visibility of this key object. We will briefly describe the 
main results and conclusions of this study which will serve as astrophysical input for the VLTI computer simulations 
presented in the next section. 

Fig. shows the SEoiliiiii and the reconstructed 2.11 ^m visibility functionEl of IRC +10420. The visibility 
0.6 at frequencies > 4 cycles/arcsec indicates, for instance, that the stellar contribution to the total flux is ~ 60% 
and the dust shell contribution is '^40%. 

An extensive grid of radiative transfer models was calculated for the dust shell of IRC +10420 assuming spherical 
symmetry and considering black bodies and model atmospheres as central sources of radiation, different silicates, 
grain-size and density distributions, various dust temperatures at the shell's inn£r boundary (determining the radius 
of the shell's inner boundary) and optical depths. We refer to Blocker et al.Ea for a full description of the model 
grid. It turned out that the observed dust-shell properties cannot be matched by single-shell models but require 
multiple components with different density distributions. The best model was found for a dust shell with a dust 
temperature of Ti — 1000 K at its inner radius of ri = 69 i?* (i?*: stellar radius). At a distance of rgw = 308 i?* 
{Y = Tsw/ri ~ 4.5) the density was enhanced by a factor of S" = 40 and its slope within the shell changed from l/r^ 
to 1/r^-^. The corresponding fits for SED and 2.11 /xm visibility are shown in Fig. |l|. The shell's model intensity 
distribution is shown in Fig. |l| and was found to be ring-like due to a limb-brightened dust-condensation zone. The 
ring diameter is equal to the inner diameter of the hot shell (^ 69 mas), and the diameter of the central star amounts 
to ~ 1 mas. 

This two-component model can be interpreted in terms of a termination of an enhanced mass-loss phase roughly 
90 yr ago. The assumption that IRC +10 420 had passed through a superwind phase in its recent history is in line 
with its evolutionary status of an object in transition from the Red-Supergiant to the Wolf-Rayet phase. 

3. INTERFEROMETRY WITH THE VLTI AND THE AMBER INSTRUMENT 

In the previous section, the spatial intensity profile of the dusty supergiant IRC +10 420 (see Fig. |^) was derived by 
means of radiative transfer models and their comparison with photometric and interferometric observations. This 
2.11 /im intensity profile will serve as object intensity profile in the simulation of monochromatic VLTI observations. 
The next steps are the simulation of light propagation from the object to the detector (through atmosphere, telescopes, 
and the AMBER wide-field mode instrument), simulation of photon noise and detector read-out noise, and finally 
data processing of the interferograms. 

3.1. Computer simulation of interferometric imaging 

Fig. ^ shows a flow chart of our simulation of interferometric imaging with the VLTI (ATs, or UTs with adaptive 
optics) and the AMBER camera in the wide-field mode (i.e. without fiber optics spatial filtering). The total 
number of detectable photons iVtot is given by the object brightness, the wavelength and band-width, the^Uecting 
area of the telescope, the total transmission of the optics, and the quantum efhciency of the detector.EEl Since 
we assume that the observations are carried out with a spectral resolution of i? = A/ A A = 70, the number of 
photons available in one spectral channel for observations in the K-hand is given by -/Vchannoi = N^.tat • Rk/R, with 
Rk = 2.2 /zm/0.4 /im = 5.5 being the spectral resolution of the A'-band (i.e. R/Rk ~ 13 is the number of spectral 
channels within the X-band; A^^^tot is the total number of X-band photons). The simulations described below refer 
to observations in one spectral channel (see also Table H) . The signal-to- noise ratio of the reconstruction can be 
increased by averaging the reconstructions of all R/Rk channels. 
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Figure 1. SED (top) and 2.11 /im visibility (middle) for superwind models with amplitudes S (=factors of density 
enhancement) at Y — rsw/ri — 4.5, as well as normalized intensity vs. angular displacement (bottom) for the 
best model with Y — r^w/ri — 4.5 and S = 40. In the bottom panel the (unresolved) central peak belongs to the 
central star. The inner hot rim of the circumstellar shell has a radius of 35 mas, and a cool component is located, 
at 155 mas. Model parameters pare: black body, T^s = 7000 K, Ti = 1000 K, tq. 55^,11 = 7.0, silicate dust grains,EZl 
standard grain size distributionEa with Omax = 0.45 /im, outer dust-shell boundary = 10** i?, (i?,: stellar radius), The 
inner shell region (Y < 4.5) obeys a 1/r^ density distribution, the outer shell region a 1/r^-^ density distribution. 
The symbols refer to the observations (see text) corrected for interstellar extinction of Ay = 5mag. [from Ref. 18] 
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Figure 2. Flow chart of the simulation of VLTI AMBER interferograms (wide-field mode, i.e. without fiber optics 
spatial filtering). 



Table 1. Parameters of VLTI- AT/ AMBER simulations: Wavelength, object brightness, optical throughput 
(^throughput of VLTI- AT + AMBER in the wide-field mode, AT mirror size: 1.8m), read-out noise and quan- 
tum efficiency of the detector, exposure time, and photon number iVchannoi per frame and per one spectral channel 
for the lowest-spectral resolution of AMBER of i? = A/AA = 70 (instrumental AMBER characteristics taken from 
Malbet et alE^I) 



wavelength 


K brightness 


opt. throughput 


read-out noise 


quant. efficiency 


exposure 


photon number iVchannoi 


2.11 


0.7 - 10.2 mag 


0.108 


15 e" 


0.6 


50 ms 


4.42 • 10^ - 7.00 • 10^ 
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In a first step, an array of Gaussian distributed random numbers is-generated and convolved with the correlation 
function of the atmospheric refraction index variations (see Roddiera) in order to generate wavefronts degraded 
by atmospheric turbulence. The typical size of the atmospheric turbulence cells is given by the Fried parameter 
ro- After the simulation of the entrance pupil, the next step incorporates the tip-tilt correction of the wavefront 
over each subpupil, but we allow for a residual tip-tilt error Stt- In the next step a typical Michelson output pupil 
is simulated by shifting both subpupils from the entrance pupil position, with a subpupil separation equal to the 
baseline, to the output pupil position where both subpupils are only separated by ~ one subpupil diameter (pupil 
reconfiguration). The output pupil is chosen such that (i) in the optical transfer function the off- axis peaks are 
separated from the central peak, and (ii) the intcrferograms are sampled with the smallest number of pixels to assure 
the lowest influence of detector noise. The following step includes light propagation through the beam combiner lens 
to the focal plane. The squared modulus of the Fourier transform of the complex amplitude in front of the beam 
combiner lens yields the intensity distribution of a Michelson interferogram of a point source. In the next step the 
required object intensity distribution is simulated (given here by the intensity distribution of IRC -1-10420) to obtain 
the Michelson interferogram of the object: The Fourier transformation of the object intensity distribution, calculated 
at those spatial frequencies covered by the the simulated interferomcter_tmseline vector, is multiplied with the off-axis 
peaks of the transfer function of the generated Michelson interferogram.EJ Finally, Poisson photon-noise and detector 
read-out noise is injected to the interferograms. The noise level depends, among other parameters, on the number 
of detectable photons, the total optical throughput p£ the interferometer (VLTI -I- AMBER in the wide-field mode) 
and the quantum efficiency of the detector. DetailsEj of the simulation parameters specific to VLTI- AT/ AMBER, 
as optical throughput, detector read-out noise and quantum efficiency, are given in Table |l|. 

3.2. Computer simulations of visibility measurements 

We performed simulations of IRC -1-10420 visibility observations in the VLTI-AT/AMBER wide-field mode and 
studied the influence of various observational parameters on the visibility accuracy. Visibility error bars were obtained 
for the following observational parameters: 

1. different seeing during the observation of object and reference star (Fried parameters ro, object and ro,rcf. ranging 
between 0.9m and 1.2m), 

2. different residual tip-tilt error (5tt, object and 5tt,rof. ranging between 0.1% and 20% of the Airy disk diameter), 

3. different object brightness (i^object^O.Tmag to 10.2mag, ifrcf.^O.Tmag). 

All computer experiments are based on 2400 interferograms of the S* = 40 intensity profile (see Fig. |l|) of 
IRC -1-10420 and refer to observations in one of the R/Rk ~ 13 spectral channels. The error bars correspond to 
the reduction of six statistically independent data sets and refer to the standard deviation a. The dependence of 
the error bar on the number of data sets has been verified (see below). In the computer experiments (1) and (2) 
mentioned above, object and reference star were assumed to have the same brightness. We chose /iT^O.Tmag (see 
Table 0) in order to minimize brightness effects. Experiment (3) simulates IRC -|-10 420 itself {K—3.5 mag) and fainter 
objects, assuming object and reference star to have the same seeing and residual tip-tilt error (ro.objcct=ro,n-fj=0.9 m; 
<5tt,object=^tt,ref=0.1%). Thcsc simulations supersede the corresponding results of a first exploratory studyE3 due to 
improved numerics and a broader statistical basis. 

The first computer experiment (different seeing conditions) consists of the simulations A to F (see top panel of 
Fig. H). Simulation A represents the case of typical seeing conditions (ro,object=0.9 m, ro,ref.=0.9 m) and an almost 
perfect tip-tilt correction (residual tip-tilt error Stt = 0.1%). It will serve as a reference for the following simulations. 
The brightness of object and reference star was chosen to be i4r=0.7mag. The visibility error ay amounts only to 
±0.0067. The remaining statistical uncertainty due to the limited number of data sets has been checked by a test 
calculation taking into account twice as many statistically independent data sets, i.e. 12 data sets instead of 6. 
The visibility error changed by 24% from ±0.0067 to ±0.0051. However, the general dependencies of the visibility 
accuracy on observational parameters (seeing, residual tip-tilt error, brightness) are only scarcely affected. 

If one considers a better seeing for both object and reference star (A-C), cry decreases in almost linear proportion 
with increasing diameter of the atmospheric turbulence cells (Fig However, improved but different seeing condi- 
tions (D-F) lead to larger visibility errors. This emphasizes that equal seeing conditions for object and reference star 
are more crucial than excellent seeing conditions of, e.g., the object only in order to obtain accurate visibilities. 
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The second set of simulations (A, G-J) illustrates the influence of a larger residual object tip-tilt error (see middle 
panel of Fig. ^). Increasing Su from 0.1% to 20% for both object and reference star increases the visibility error from 
0.0067 to 0.0084 (A,F,G). Like in the simulations of different seeing conditions, it is essential to have the same tip-tilt 
error for both object and reference star in order to achieve high accuracy. Differences in Stt enhance the visibility 
error (H,I). 

Finally, in the third experiment of this series an IRC -1-10 420-like intensity distribution is assumed but much 
much fainter (i.e. fainter than iir=3.5mag) objects are considered (see bottom panel of Fig. ^, simulations A, K-N). 
The visibility error stays almost constant up to an K magnitude of 7 mag. If the simulated object if-magnitude is 
9.2 mag, the visibility error has increased by a factor of ~ 2, i.e. to ay = 0.0156, but is still acceptable although the 
photon number decreased to A^channci = 1-76 • 10'^ (i.e. by a factor of 2500, see Table |^). Fig. || demonstrates that 
for even fainter objects noise becomes more and more important leading to an steep increase of the visibility error. 
However, the signal-to-noise ratio can be improved (~ R/ Rk ^ 3.6) if the visibility measurements of all R/Rk 
spectral channels within the if -band are averaged. 

For illustration, Fig. ^ shows the results of the simulations A*, G*, J* and D* for baselines of 50 and 100m 
together with the model predictions for IRC -1-10420. The asterisk indicates that within these simulations the object 
brightness is that of IRC -f 10420, viz. 1^=3.5 mag (photon number: 3.47 ■ 10^) As it is obvious from Fig. ^ (bottom 
panel), this brightness is well within the range where the visibility error is not very sensitive to brightness changes. 
Accordingly, the results are very close to those of the simulations A, G, J, and D. Again, different seeing conditions 
(simulation D*) for object and reference star are more crucial than, e.g., residual tip-tilt errors (simulations G* and 
J*). For typical seeing conditions (Fried parameter differences ^ 10%), the visibility accuracy achievable with the 
VLTI using the ATs in the wide-field mode is clearly sufficient to distinguish between different radiative transfer 
models of IRC -f 10 420 and, thus, to prove (or disprove) theoretical predictions. 

3.3. Computer simulations of phase-closure imaging 

In addition to the above visibility studies, image reconstruction simulations were performed as well. For this purpose, 
modulus (—visibility) and phase of_the Fourier transform of the object's intensity distribution were determined 
utilizing the phase-closure method.u The phase closure method requires the simultaneous measurement of the 
object with at least three telescopes. For a corresponding computer experiment we chose three telescopes covering 
five configurations. The baseline between telescope 1 and telescope 2 was 8 m whereas the baseline between telescope 
2 and telescope 3 amounted to 8, 16, 24, 32 and 40 m in order to facilitate a simple recursive algorithm for phase 
reconstruction from the measured closure phases. A non-redundant arrangement of the output pupil ensures a 
sufficient separation of the object information in frequency space. 

As model object served a binary star with an intensity distribution as shown in Fig. |^ (bottom; dashed line). 
The components' intensity ratio is 1:2. As in the case of the visibility simulations of IRC -1-10 420 each simulation is 
based on 6 statistically independent data sets. The number of interferograms per data set is 1200. 

Fig. ^ shows the simulated modulus, phase and image reconstruction of the binary for total object /^-magnitudes 
of 0.7 mag and 9.2 mag, resp. Even for objects as faint as K = 9.2 mag the agreement between test object and 
reconstructed object is very good. Beyond K = 10 mag the errors increase considerably. This is also illustrated in 
Fig. ^ which shows the photometry error (i.e. the deviation from the components' intensity ratio of 1:2) as a function 
of the object brightness. The signal-to-noise ratio of the reconstruction can be increased if the reconstructions of all 
R/Rk spectral channels are averaged. 

4. CONCLUSIONS 

We have presented computer simulations of interfcrometric imaging with the VLT interferometer and the AMBER 
instrument in the wide-field mode. These simulations include both the astrophysical modelling of a stellar object 
by radiative transfer calculations and the simulation of light propagation from the object to the detector and sim- 
ulation of photon noise and detector read-out noise. We focussed on stars in late stages of stellar evolution and 
examplarily studied one of its most outstanding representatives, the dusty supergiant IRC -1-10 420. The model in- 
tensity distribution of this key object, obtained by radiative transfer calculations, served as astrophysical input for 
the VLTI/ AMBER simulations. The results of these simulations show the dependence of the visibility error bar on 
various observational parameters. With these simulations at hand one can immediately see under which conditions 
the visibility data quality would allow us to discriminate between different model assumptions (e.g. the size of the 
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Figure 3. Dependence of the error bars of simulated AT-VLTI/AMBER (wide-field mode) observations of IRC +10 420 at 
2.11 /im on the seeing conditions [top], residual tip-tilt error [middle], and object brightness [bottom]. The table columns refer 
to the Fried paramater ro, the residual tip-tilt error Su, the AT- magnitude and the visibility error cry (based on 6 statistically 
independent repititions of each simulation). Each AT-VLTI/AMBER simulation refer to A=2400 interferograms of the S=40 
intensity profile (see Fig. |l|). j 
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Figure 4. Dependence of the error bars of simulated AT- VLTI/ AMBER (wide-field mode) observations of 
IRC +10420 at 2.11 /j,m on (i) seeing differences of object and reference star observations (Fried parameter, ro, 
differences; simulation D), and (ii) different residual tip-tilt errors, ^tt, of object and reference star observations (sim- 
ulations G and J). Lines refer to radiative transfer models of different superwind amplitudes S. The symbols refer 
to AT- VLTI/ AMBER simulations (each with 7V=2400 interferograms) of the 5=40 intensity profile (see Fig. |l|). for 
baselines of 50 and 100 m (100 m: q—227 arcsec^^). To better distinguish between the simulations, the data points 
belonging to one baseline are somewhat shifted with respect to the spatial frequency. The error bars are based on 6 
statistically independent repetitions of each simulation. The table gives the parameters of the simulations A*, G*, J*, 
and D* (cf. Fig. ^). The asterisks indicate that the object brightness is here 3.5 mag as appropriate for IRC -1-10 420 
whereas the simulations shown in Fig. 13 refer to an object brighness of 0.7 mag. 
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Figure 5. Modulus (top), phase (midlle) and image (bottom) of simulated AT- VLTI/ AMBER (wide- field mode) phase- 
closure observations of a binary at 2.11 ^m if three telescopes are employed (reconstructions are derived from interferograms 
of only one of the R/Rk ~ 13 spectral channels). The left column refers to a total object brightness oi K = 0.7mag, the right 
column to one oi K = 9.2 mag. The AT- VLTI/ AMBER simulation takes 5 different telescope configurations into account and 
is based on 6 statistically independent data sets, each consisting of 1200 interferograms. The crosses and solid lines refer to 
reconstructions based on one data set, the dashed lines to the original object. 
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Figure 6. Dependence of tlie error bars of simulated AT- VLTI/ AMBER (wide-field mode) phase-closure observations of a 
binary at 2.11 fim on the total object brightness if three telescopes are employed (reconstructions are derived from interferograms 
of only one of the R/Rk ~ 13 spectral channels). The table columns refer to the Fried paramater ro, the residual tip-tilt error 
Stt, the A'-magnitude and the photometry error (jp (i.e. the deviation from the components' intensity ratio of 1:2; based on 6 
statistically independent repititions of each simulation). Each AT- VLTI/ AMBER simulation refers to A=1200 interferograms 
of a binary with an intensity profile as given in Fig. U. 



superwind amplitude S). Inspection of Fig. |j shows that in aU studied cases the observations will give clear preference 
to one particular model. Therefore, observations with VLTI will certainly be well suited to gain deeper insight into 
the physics of dusty supergiants. 
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